Phase transitions and spectral properties of the ionic Hubbard model in one dimension 
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The ionic Hubbard model is investigated at half filling at zero temperature. We apply the cellular 
dynamical mean-field theory to the one-dimensional ionic Hubbard model to compute local quantities 
such as double occupancy and staggered charge density. Both quantities provide general transition 
behavior of the model from a band insulating phase to a Mott insulating phase. The renormalized 
band gap is introduced as an efficient order parameter for the transition from a band insulator. We 
also present the spectral properties of the ionic Hubbard model, which exhibit characteristic features 
for both weak and strong interactions. 



I. INTRODUCTION 

Strongly correlated electron systems have been one of 
the most interesting problems in modern condensed mat- 
ter physics. The interest in strongly correlated systems 
was motivated by experiments on transition-metal oxides 
which were inconsistent with the predictions of the con- 
ventional band theory. Since the argument by Mott that 
the mutual interaction between electrons can cause insu- 
lating behavior, extensive research has been performed 
on the effects of the interaction in the metal-insulator 
transitionii It is practically impossible to include all the 
degrees of freedom in strongly correlated systems, and 
the standard approach has been to solve theoretical mod- 
els constructed with essential ingredients. The Hubbard 
model (HM) is one of the most popular models in strongly 
correlated systems;^ It includes only two essential com- 
ponents: electron hoppings and local Coulomb interac- 
tions. Although the model seems extremely simplified, it 
has successfully described the metal-insulator transition 
caused by the mutual interaction between electrons. 

The ionic Hubbard model (IHM) , an extended version 
of the HM, was proposed to explain the neutral-ionic 
transition in the quasi-one-dimensional charge-transfer 
organic materials)^"— Unlike the original HM, this model 
is an insulator in the absence of the mutual Coulomb 
interaction. It enables one to examine an interesting 
insulator-insulator transition, from a band insulator (BI) 
under weak interaction to a Mott insulator (MI) under 
strong interaction. Accordingly the transition nature in 
the IHM is expected to differ significantly from the usual 
metal-insulator transition in the HM. 

A more interesting feature of this model is the possi- 
bility of a nontrivial intermediate state, sandwiched be- 
tween the two insulating phases. Extensive studies of the 
intermediate state in the IHM have been carried out in 
various spatial dimensions. In the infinite dimensions, 
which can be treated exactly by the dynamical mean- 
field theory (DMFT),— the metallic phase is observed as 
an intermediate phase for the weak alternating potential, 
whereas a direct insulator-insulator transition is shown 



under a strong staggered potential.—"— The same con- 
clusion has been obtained by a recent study using a co- 
herent potential approximation^^ In two dimensions, on 
the other hand, there has been some controversy as to the 
nature of the intermediate phase in the IHM. Quantum 
Monte Carlo calculation a^"^!^^ showed that the metallic 
phase exists, as is observed in the infinite dimensions. 
In contrast, the cellular dynamical mean- field theory 
(CDMFT)i^ as well as the variational cluster approacbi^ 
predicts a bond-ordered insulating phase as an interme- 
diate state. 

In one dimension, which is our main interest in this 
paper, the possibility of an insulating intermediate state 
was suggested by the bosonization methodfii which pre- 
dicts a spontaneously dimerized insulating (SDI) phase 
between BI and MI phases. Many subsequent inter- 
esting works for the one-dimensional (ID) IHM have 
been reportedJ^"— Particularly density matrix renor- 
malization group (DMRG) calculations have confirmed 
the existence of an SDI phase for intermediate interac- 
tion strengths jiSii^ However, the critical properties of 
the transition between BI and SDI phases obtained by 
DMRG calculations turned out not to be consistent with 
those of two-dimensional (2D) Ising universality class 
predicted in a bosonization method. For the transition 
between SDI and MI phases, even the critical interaction 
strength has not been identified clearly from finite-size 
scalings. 

The main purpose of this paper is to investigate in de- 
tail the nature of transitions in the ID IHM at half filling. 
In order to achieve this we use the CDMFT)^ which is 
one of the cluster extensions'^ of DMFT. While a sin- 
gle site is chosen to construct the self-consistent equa- 
tion in the DMFT, the CDMFT picks up a cluster com- 
posed of several sites in the spatial dimension of the sys- 
tem. This makes it possible to include short-ranged spa- 
tial fiuctuations inside the cluster, which are expected to 
be important in low-dimensional systems. The CDMFT 
turns out to be a very efficient method even in one di- 
mension,—"— the worst case for a mean-field theory. In 
particular, it has been shown that the CDMFT of the 
ID HM yields excellent agreement with the Bethe ansatz 
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exact solution<21^ Such a good agreement in turn re- 
inforces our expectation of an accurate description by 
the CDMFT method of the ID IHM for which no exact 
solution is available. We employ the exact diagonaliza- 
tion method as an impurity solver to study the ground 
state of the ID IHM. The exact diagonalization method 
is powerful since it deals with all the quantum fluctua- 
tions on an equal footing, although it limits the number 
of sites inside the cluster as well as that of bath sites. Re- 
cent studies on correlated systems show that single-site or 
cluster DMFT combined with the exact diagonalization 
gives reliable results for finite temperatures as welli^"— 

We calculate local quantities such as staggered charge 
density and double occupancy, which have been reported 
to be evaluated accurately in the CDMFT approach<22i^ 
We also introduce a renormalized band gap defined by 
the self-energy corrected band gap and demonstrate that 
it plays the role of convenient order parameter for the BI 
phase. Any indication of the transition from an interme- 
diate phase to an MI phase is not observed in the local 
quantities, such as the abrupt change in slope which was 
reported in the cellular dynamical mean-field study of the 
two-dimensional IHM.— In order to gain further insight 
into the transition nature we finally focus on the spectral 
properties of the model. Computed spectral weights re- 
veal the spin-charge separation in an MI phase, which is 
characteristic of the ID system. Such spin-charge separa- 
tion has been reported by earlier studies through cluster 
extensions of the DMFT only for the ID EM^^ The 
transition from an intermediate phase to an MI phase 
exhibits rather a crossoverlike behavior, which may be a 
reason for the difficulty in obtaining the transition point 
clearly in earlier numerical investigations. 

This paper is organized as follows. In Sec. [Hi we de- 
scribe the IHM and introduce briefly the procedure of 
the CDMFT. The numerical results and discussions are 
presented in Sec. HID We give a summary in Sec. IIVI 



II. MODEL AND METHOD 



The Hamiltonian of the ID IHM is given by 
H = -tY^ (c.Lsa + H.c.) ""^ 

{i.j)cr la 



(1) 



where (cj^) creates (destroys) an electron with spin 
a at the tth site and n^^=c^^^c^^. The hopping of elec- 
trons is allowed only between the nearest neighbors. The 
parameters i, J7, and ^ are the hopping amplitude, the 
on-site Coulomb repulsion, and the chemical potential, 
respectively. Throughout the paper we will represent all 
the energies in units of t. The system is composed of two 
alternating sublattices A and B. The potential energy 
difference between nearest neighbors is A(> 0) and every 



site belonging to the sublattice A (B) has lower (higher) 
potential energy by A/2 than the chemical potential. If 
A vanishes this model is restored to the original HM. We 
set the chemical potential to be half the Coulomb repul- 
sion {fi ~ U/2) in order to maintain the half filling of 
electrons over the system. 

We use the CDMFT to study the ID IHM. The infinite 
lattice is reduced to a cluster of size Nc which hybridizes 
with the self-consistent electronic bath sites within the 
CDMFT. We employ the exact diagonalization method 
as an impurity solver. In order to obtain the cluster 
self-energy, we map the lattice model to the impurity 
Hamiltonian, 

^y<7 fj, 

IJ.la la 

(2) 

where /i, = 1, 2, . . . , Nc are the cluster site indices and 
I = l,2,...,Ni, label the bath sites. The matrix E 
contains the hoppings, the staggered potential, and the 
chemical potential inside the cluster, which is given ex- 
plicitly by 



E : 



/-/i - A/2 -t 

-t -//-A/2 



■•■ -/i-f A/2 / 

(3) 

[The circumflex over a symbol (") represents a. Nc x Nc 
matrix.] The parameters {eia} and {Vf^ia}, which re- 
spectively denote the bath energy levels and hybridiza- 
tion strengths with the clusters, are determined from the 
imposed self-consistency conditions. From the impurity 
Hamiltonian we compute the cluster Green function G as 
well as the cluster self-energy, E'^ = — G~^, where Q 
is the Weiss fleld describing the noninteracting bath. 

The Green function of momentum k in reduced Bril- 
louin zone is given by 



G{k, iujn) = M - i{k) - S'=(ia;„) 



(4) 



where M ^j^y=[iuin — E^^)S^i,, i{k) is the Fourier trans- 
form of the hopping matrix, and a;„ = {2n + 1)tt//3, 
n = 1,2,..., iVniax are fictitious Matsubara frequencies. 
Here we used /3 = 100 and iVmax = 400 and omitted 
the spin index a for simplicity. The local lattice Green 
function Gioc is then determined by 



Gloc(iW„) 

The new Weiss field Gnaw 
consistent equation. 



k 



G{k,iujn). 



(5) 



is obtained through the sclf- 



(6) 
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the best fit of which produces new bath parameters 
{ezcr, V^io-}- The above procedure is repeated until the 
convergence is reached. For a more detailed procedure, 
we refer the reader to earlier works»22iM 

It is known that even the cluster of iVc = 2, which is the 
minimal size describing the system within the CDMFT 
approach, is sufficient to obtain an accurate estimation of 
the local quantities in the ID HM^^i Further, we have ob- 
served that in the ID IHM the local quantities obtained 
by the calculations for various combination of Nc and Nb 
are not significantly dependent on the choice of Nc and 
A^h, although larger Nc has some tendency to improve the 
momentum resolution of spectral weights. The number 
of sites in the impurity Hamiltonian is practically limited 
hy Nc + Nb< 12. We use iVc = 4 and TV^ = 8 for the pre- 
sentations of most quantities investigated in this paper, 
and other choices of Nc and N^ will be specified. 



III. RESULTS AND DISCUSSION 

A. Local quantities 

In order to examine how the system evolves with the 
variation of the local interaction we first calculate two 
local quantities: staggered charge density and double oc- 
cupancy. The staggered charge density is defined by the 
difference between the number densities at two sublat- 
tices, nA~'nB^ where the sublattice number densities can 
be calculated as 
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(7) 



with the angular brackets being the expectation value 
over the ground state of the impurity Hamiltonian. We 
also calculate the double occupancy defined by 



(8) 



which is known to be a convenient measure in the tran- 
sition from metal to an MI phase. 

In the ionic limit (A ^ t,U)^ it is energetically favor- 
able that all the electrons are in the sublattice A, pro- 
ducing unity of the staggered charge density. As U is 
increased, the energy cost of two electrons to stay in the 
same site becomes large, reducing the staggered charge. 
In the strong-coupling limit {U » A), the staggered 
charge is expected to approach zero. 

In the case that the system has neither the interaction 
nor the staggered potential {U = A = 0), the double 
occupancy is 1/4 since every site has an equal possibility 
to be occupied independent of spin. The presence of the 
staggered potential in the system tends to increase the 
double occupancy. In contrast, the repulsion between 
electrons makes the doubly occupied sites less favorable 
and results in the reduction of the double occupancy. 





FIG. 1: (Color online) (a) Staggered charge density nA—ns- 
(b) Double occupancy D as a function of U/t for A/t = 1. 
Both quantities approach zero monotonically as U is in- 
creased, signifying that the system exhibits a BI phase for 
weak interaction while an MI is recovered in the limit of strong 
interaction. The lines are merely guides to the eyes. 



The features of two quantities described in the above 
reasoning are generally consistent with our numerical re- 
sults in Fig. [1] which shows (a) the staggered charge 
density and (b) the double occupancy as a function of 
U . The staggered charge is relatively large in the weak 
interaction regime and monotonically decreases with in- 
creasing U ; it reveals that the system exhibits a BI phase 
for weak interactions and that an MI behavior turns up in 
the regime of strong interactions. In our numerical work 
we have also confirmed that the system with large A 
shows larger staggered charge density over the whole re- 
gion although the increment is reduced for stronger inter- 
actions (not shown). The analysis of double occupancy 
has drawn a similar conclusion. In the noninteracting 
system {U = 0) the double occupancy is larger than the 
free-electron value 1 /4 and the system lies in a BI phase. 
The double occupancy is reduced monotonically by the 
increase of U , approaching zero, as is expected in an MI. 
Monotonic decrease in both the staggered charge density 
and the double occupancy demonstrates that the increase 
of the interaction strength drives the system from a BI 
into an MI. It is of interest to note that no abrupt change 
is observed in the variation of both quantities with an in- 
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FIG. 2: (Color online) Energy densities A/t — 1 as a func- 
tion of U. The total, kinetic, potential, and correlation energy 
densities are denoted by squares(n), circles(O), triangies(A), 
and diamonds(O), respectively. The corresponding energy 
densities in the ID HM are given as lines for comparison ex- 
cept for the potential one which is zero in the HM. 



crease of the interaction strength. 



B. Energy densities 

In order to gain more insight on the transition between 
BI and MI, we compute the energy density of the system. 
The ground-state energy density of the system is calcu- 
lated as 



E 
N 



ex + CA + ec/, 



|EEf]^Tr{£(fc)G(fc,*.;„)} 



(9) 



where e/^, a-nd eu denote kinetic, potential, and cor- 
relation energy densities, respectively. The factor 2 in 
€if comes from spin degeneracy. Computed energy den- 
sities illustrate which contribution plays a dominant role 
in each phase. In Fig. [2] we have plotted various con- 
tributions to the energy density of the ID IHM along 
with those of the ID HM for comparison. Since the HM 
does not have staggered potential, the contribution from 
potential energy is zero and is not plotted in the figure. 

In the region of weak interactions {U < 2t), we can see 
that both the kinetic and the correlation energy densities 
are higher than those of the standard HM. The energy 
gain in the potential contribution compensates for the 
increase in other energy densities, demonstrating the BI 
nature of the system in this region. The increase of the 
interaction strength lessens such a tendency gradually. 
In the strong-interaction region {U > At) all the energy- 
density contributions are almost the same as those of the 



0.4 

0.2 


0.4 

0.2 


0.4 

0.2 


0.4 

0.2 


0.4 

0.2 





(a) 



(b) 




(c) 



(d) 



(e) 



-10 





(H/t 



10 



FIG. 3: (Color online) Local density of states p{uo) for 
A/t=1.0 and U/t^O.3, 2.4, 6.0, 8.0, and 10.0 from top to bot- 
tom. In the weak-coupling regime the system displays a two- 
band structure and the gap between the two bands is reduced 
gradually with the increase of U. In the strong-coupling limit, 
on the other hand, the two bands move away from the Fermi 
level and each band seems to be split into two subbands. A 
broadening factor £=0.05 is used. 



HM and the energy contribution from the staggered po- 
tential is almost zero, signifying that the system is in an 
MI phase. As in the local quantities examined in the pre- 
vious subsections, all the energy densities display contin- 
uous variations with the increase of interaction strength. 



C. Local density of states 

The local density of states (LDOS) provides more 
detailed information on the single-particle properties. 
Within the CDMFT approach, the LDOS is given by 



(10) 



where A{k,uj) is a spectral weight with the energy uj and 
the momentum vector k in the full Brillouin zone. We re- 
store the translational symmetry broken in the CDMFT 
formalism by the periodization of the Green function for 
each sublattice 



Gaik,LL>) 



-ik{}i — i>) 



G(fc,w)]^, (a = A,B), 
(11) 
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where G{k,u;) is given in Eq. ([4]). Then we can compute 
the total spectral weight yl(fc,aj) as 

A{k,uj) = -- V lmGa{k,uj + ie), (12) 

TT ^ — ' 

a=A,B 

where Ga is the periodized Green function of the sublat- 
tice a and e is a small broadening factor. In this work 
we used e=0.05. 

The LDOS is shown in Fig. [3] for several values of U. 
In the noninteracting system {U = 0), the LDOS can 
be computed analytically and is composed of two bands 
which are separated by a band gap A due to the stag- 
gered potential. The CDMFT results generally repro- 
duce the analytical LDOS for a noninteracting system as 
demonstrated in Ref. [H. Turning on the interaction does 
not change the LDOS very much from the noninteract- 
ing LDOS, and the two band structure is retained. For 
weak interactions the increase of interaction strength re- 
duces monotonically the band gap around a Fermi level. 
Around a certain value [/o~2.4t the LDOS displays a 
minimum spectral gap and prominent long tails show up 
at the outer edges of the bands. Further increase of U 
above Uq in turn enlarges the gap between the two bands. 
Each band is apparently split into two subbands and the 
LDOS is composed of four bands for strong interactions. 
Such a four-band structure is also observed in the original 
HM, where it is caused by spectral weights concentrated 
on the spinon-holon continuumi^ In the IHM the gap 
between the subbands is distinguished more clearly than 
in the HM and it is found to be proportional to A. We 
will give more detailed discussions on this topic in the 
section that deals with the spectral weights. 

D. Spectral gap and renormalized band gap 

For a quantitative analysis of the gap around a Fermi 
level, we investigate the spectral gap d, which is defined 
as the energy difference between the highest filled and the 
lowest empty levels in the LDOS. Since the IHM always 
displays a minimum gap at k—T:/2 the spectral gap can 
be conveniently obtained from the spectral weight at the 
Fermi point k~'K/2 within the CDMFT. Since the sys- 
tem has a particle-hole symmetry, we estimate the spec- 
tral gap 5 to be twice the difference between the Fermi 
level and the peak of w) closest to it for positive w. 
The measured gap is plotted as a function of the inter- 
action strength U for various A in Fig. [H The overall 
behaviors of the spectral gap are consistent with those 
observed in the LDOS. The spectral gap is equal to A in 
the noninteracting system (C/=0). It is reduced from the 
noninteracting value A by weak interactions, reaches a 
minimum at [/ = Uq, and increases with U in the regime 
of strong interactions above Uq. Our CDMFT results of 
the system with various A in Fig. |4]have shown that Uq 
increases with the bare band gap A. 

In a BI quasiparticle excitations are well defined and 
we can interpret the spectral gap as the energy difference 
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FIG. 4; (Color online) Spectral gaps 5 as a function of the 
interaction strength U for various A/t. The spectral gaps 
for the system with A/t = 1,2,3, and 4 are denoted by 
diamonds(O), squares(n), circles(O), and triangles(A), re- 
spectively. With the increase of U the spectral gap decreases 
for weak interactions while it grows larger in the regime of 
strong interactions. 
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FIG. 5: (Color online) Difference between the spectral gap 
5 and the renormalized band gap Aren for various values of 
A/t. The same symbols are used as in Fig. [l] 



between the lowest quasielectron and quasihole excita- 
tions. By expanding the self-energy for each sublattice 
around the Fermi level we can estimate the position of 
the poles of the Green function on the real frequency 
axis near the Fermi level. We then define the renormal- 
ized band gap Aron by the difference between the poles 
of the Green's function of the sublattice A and B which 
are the closest to the Fermi level, which is expected to 
be the same as the spectral gap in a BI phase. By defini- 
tion Arcn equals a bare band gap A in the noninteracting 
system since the self-energy vanishes and the poles of the 
Green function are identical to the bare dispersion of the 
system. Indeed Aj-en is a band gap with a self-energy cor- 
rection. On the other hand, in correlated phases such as 
an MI, the quasiparticle excitations are not well defined 
and the discrepancy between 5 and Aion will show up. 
We can thus expect that the comparison of the renor- 
malized band gap Arcn with the spectral gap 5 will be a 
good parameter for determining whether the system is in 
a BI phase or not. 
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TABLE I: Comparison of the critical interaction strength Ud 
below which the system is in a BI phase and Uo at which the 
system shows a minimum spectral gap. All the energies are 
given in unit of t. 



A 


Uo 


f7ci 


1.0 


2.4 


2.2 


2.0 


3.4 


3.2 


3.0 


4.4 


4.1 


4.0 


5.3 


5.0 



Within the CDMFT we have computed and plotted 
the difference between the spectral gap 6 and Aj-cn in 
Fig. [S] We can see that it is zero for C/=0, which is 
guaranteed from the definition of the renormalized band 
gap. It is remarkable that the difference remains zero 
over a finite region of U below a certain critical value 
Uci- Above Uci the discrepancy between S and A.^on 
turns out to grow rapidly, signifying that an MI phase 
or another correlated phase emerges for U > Ud- For 
A = t, the CDMFT yields Ud « 2.2t; this is compara- 
ble to, although slightly smaller than, the critical values 
obtained in the existing works, Ud ~ 2.67i from the 
DMRG studyi^ and Ud ~ 2.3i from the effective model 
in the strong-coupling limiti^ We compare Ud and Uq 
for various A in Table ID It is found that two interaction 
strengths correlate very much with each other and Uq is 
always slightly higher than Ud-, which is also consistent 
with existing worksi ^^'^^ 



E. Spectral weights 

A useful quantity for demonstrating momentum- 
resolved information on the correlation effects in the sys- 
tem is a total spectral weight A{k,Lu) given in Eq. (fT^ . 
This is computed from periodizcd sublattice Green func- 
tions in Eq. (fTT]) and we plot A{k,uj) for various interac- 
tion strength U in Fig. [6l 

For weak interaction below Ud two quasiparticle bands 
are sharply defined in A{k,uj) as illustrated in Fig. [6l^a) 
for U=0.3t. Such quasiparticle bands which are sepa- 
rated by a band gap result in two-band structure of the 
LDOS which is similar to that of a nonintcracting system. 
In this regime only the gap between the two bands and 
the band widths are renormalized by weak interactions. 

In Fig.m^b) we can observe that around Ud the single- 
particle dispersion begins to be broadened particularly 
around the zone center fc=0, which is a signature of 
the transition to a correlated phase from a BI phase. 
The broadening around the zone center is a source for 
the appearance of long tails at the outer edges of the 
bands in the LDOS. Here it is of interest to note that 
the broadening shape of the dispersion closely resembles 
the spin-charge separation observed in the ID HM. The 
spinon (holon) is a fractionalized excitation which car- 
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FIG. 6: (Color online) Spectral weights for A/t — 1.0 and 
(a) U/t = 0.3; (b) U/t = 2.2; (c) U/t = 6; and (d) U/t = 8. 
We used a cluster of size Nc=6 with Nb—G bath sites. Since 
the unit cell of the IHM is twice that of the original HM, the 
first Brillouin zone is reduced to — 7r/2 < k < 71/2. For better 
resolution the spectral weights are rescaled according to their 
maximum values in each plot. 



ries only spin (charge) but no charge (spin). In the ID 
HM, the Bethe ansatz solution gives the dispersion of the 
exotic excitations. The assembled dispersions construct 
the spectral weight and we can observe some prominent 
dispersions such as the spinon and the holon branches. 
In the ID HM most weights of A{k, ut) are concentrated 
on the spinon branch under relatively weak interaction 
and more weights transfer to the holon branches and 
other accessible regions of higher energies as the inter- 
action strength is increased)^ Similar weight transfer is 
observed in the spectral weights of ID IHM as shown in 
Fig.El 

Although the spectral weights of the two models are 
quite similar to each other in the strong interaction 
regime, a remarkably different feature also arises from 
the presence of staggered potential. In the absence of 
the staggered potential two holon branches for fc>0 and 
for fc<0 cross each other at the zone center slightly above 
uj=U /2, and extends as secondary holon branches in the 
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FIG. 7: (Color online) Spectral weight for Nc ^ 6, Nt = 6, 
U/t = 10.0, (a) A/t = 0, (b) A/t = 2, and (c) A/t = 4. 



other region, as is reproduced in Fig.[7[a)i^ The spectral 
weights for A/t = 2 and 4, which are shown in Figs.[7l^b) 
and[7l^c), demonstrate marked gaps at the crossing points 
of the holon branches. The fact that the gap width is pro- 
portional to A also supports that the degeneracy of the 
holon branches at k=0 is lifted by the staggered potential. 
Accordingly, the dispersion displays four well-separated 
bands, yielding the characteristic four-band structure ob- 
served in the LDOS for large U in Fig. [3l We also note 
that the shift of the spinon and the first holon bands 
toward the Fermi level produces some reduction in the 
Mott gap in the presence of a staggered potential. 

The comparison of the spectral weights with different 
Nc gives us a good guide to the understanding of the 
overall distribution of spectral weights. In Fig. [8] we plot 
the spectral weights on the two momentum points k=0 
and fc=7r/2 with different cluster sizes Nc=2, 4, and 6. In 
the plot of spectral weights at the zone center in Fig. [H^a) 
we can recognize three prominent peaks around uj « 3t, 
4.3t, and 5.5t, which are identified as the spinon and the 
two split holon branches, respectively. A similar ten- 
dency is demonstrated clearly in the plot of A{t:/2, uj) in 
Fig. [Hljb). The first peak around uj k, 2t corresponds to 
a merging band of the spinon and the lower holon band, 
while the higher holon band generates the second peak 
around w « 6.5f. The peaks which seem to be prominent 
for Nc=2 tend to be smeared out to a continuum with 
the increase of Nc- 

Since the pioneering bosonization approacf^ it has 
been generally believed that the ID IHM shows two suc- 
cessive transitions with an increase of the interaction 



I — — ■ — ^ 1 

2 4 6 8 

co/f 

(b) 

FIG. 8: (Color online) Spectral weight for A/t = 1, U/t = 8, 
(a) k = Q, and (b) k = 7r/2 with different cluster sizes. 



strength. In the CDMFT we have identified the first tran- 
sition from the BI by examining the rcnormalizcd band 
gap. However, we are not able to position the second 
transition point from the intermediate insulating phase 
to the MI in the investigation of the systems with var- 
ious values of A. We have not observed any nonana- 
lytic behaviors in local quantities such as the staggered 
charges or the double occupancy, in contrast to the 2D 
IHM where some kinks in the local quantities were pro- 
posed as a signature of the second phase transitions^ No 
abrupt change in the spectral properties occurs as the 
interaction strength is increased. Particularly the spec- 
tral weights exhibit a rather gradual transition to the MI. 
Although the origin of such difficulty in positioning the 
second phase transition is not clear, we believe that it is 
related to rather slow falloff of the bond-order parameter 
around the second transition in the DMRG studyJ^ 



IV. SUMMARY 

We have investigated the one-dimensional half-filled 
ionic Hubbard model at zero temperature using the cellu- 
lar dynamical mean-field theory. We have computed the 
staggered charge density and the double occupancy. Both 
quantities display monotonic decrease with an increase 
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of the interaction strength, signifying that the system 
evolves from a band insulator to a Mott insulator. The 
energy-density analysis shows that the potential energy 
gain occurs for weak interactions while the system with 
strong interactions gives almost the same energy contri- 
butions as in the Hubbard model. The phase boundary of 
a band insulating phase has been determined by the com- 
parison of the rcnormalized band gap with the spectral 
gap. Around the phase boundary we have also observed a 
minimum of the spectral gap for various strengths of stag- 
gered potential. We have calculated the spectral weights 
and analyzed the detailed structure of spin-charge sep- 
aration by analogy with the one-dimensional Hubbard 



model. The staggered potential turns out to produce a 
gap at the crossing point of two holon branches, which is 
proportional to a band gap parameter. 
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